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Abstract 

b ■ 

0^ I We derive the bulk viscous damping timescale of hybrid stars, neutron stars with 

quark matter core. The r-mode instabihty windows of the stars show that the theo- 
. retical resuhs are consistent with the rapid rotation pulsar data, which may give an 

indication for the existence of quark matter in the interior of neutron stars. Hybrid 
PLh' stars instead of neutron or strange stars may lead to submillisecond pulsars. 
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Probing the composition of matter in the interior of compact stars is of in- 
terest for understanding of pulsars. The available theories have not yet given 
a unique solution to this problem. Since the deconfinement transition from 
hadron matter into quark matter is possible at high density, one logically pro- 
poses that some pulsars could be strange stars. But how to distinguish strange 
stars from neutron stars is a difficult and an urgent issue. First, one suggested 
that bulk properties of neutron stars such as radius and moment of inertia, 
may be marked off due to the difference between self-bound and gravity-bound 
matter. Unfortunately, the bulk properties of t hese two kinds of stars are simi; 



lar in the observed mass range 1 < M/Mq < 2 fjHaensel. Zdunik and Schaeffer 



19861 ). Another way is to look at the possible difference of their highest rotation 
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frequency, since the discovery by Andersson, Friedman and Morsink of r-mode 
instabilities in the rotating compa ct stars put r ather severe hmits on the high- 



est r otation frequency of pulsars (j Andersson I . Il998l : iFriedman and Morsink 



19981 ) ■ A pioneer investigation found that the limiting rotation of strange stars 
is in milliseco nd period range but fo r normal neutron stars in the period range 



of 15- 20ms flMadsen I . Il998l . l2000h . However, a strange star evolving to the 

limiting frequency corresponds to temperature in 10® 10''K range instead 

of the inferred core temperature of millisecond pulsars of a few times lO^K. 
Obviously, this is a serious drawback if we attributed the uncorrespondence to 
the under-r epresenting stat istics due to a low number of objects as suggested 



by Madsen ljMadsen 1 . 12000| ) . Although some works also argued that neutron 



stars with solid crusts would be the candidates to explain the rapidly rotating 
accreting pulsar's data through thermal runaway recycle of r-mode. Due to 
the p resence of viscous bou ndary layer damping (IBildsten and Ushomirsky 



tne presence oi viscous Dou naary layer dampme iimidsten ana usnomirsKy . 
2000l : lAndersson eraP . l2000h . the resulting viscous heating is so intense that 



it can heat the crust-core interface to the melting temperature of the solid 
crust when the r-mode amplitude is larger than some critical value, which is 
crudely estimated to be ~ 10~^ bvlOwenI (119991) and a naqre ac curate value is 



given as 5 X 10 ^ by (iLindblom. Owen and Ushomirsky 1 . 12000| ) 



Following the above train of thought, we will here investigate the r-mode 
instability window of hybrid stars, i.e., the neutron stars containing quark 
matter core. We show that the instability window of hybrid stars explains the 
millisecond pulsar data very well, but the strange stars and normal neutron 
stars would not fit because of the weaknesses concerned above. 

In rotating relativistic stars, gravitational radiation drives the r-mode while 
various dissipation mechanisms counteract the fluid motion. In general, shear 
viscosity, or surface rubbing that could also be applied to the hybrid stars, 
which also have a solid crust similar to the normal neutron stars, suppresses 
the mode at low temperatures but bulk viscosity dominates at high temper- 
atures. The critical rotation frequency for a given stellar model as a function 
of temperature follows 
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where Tgr < is the characteristic timescale for energy loss due to gravity wave 
emission, Tg^ and Tbv are the damping times due to shear and bulk viscosities, 
and Tgr is that for surface rubbing due to the presence of the viscous boundary 
layer. 



A polytropic equation of state with a low in dex n is a good approximat i on for 



compact stars, as discussed in many papers ijKokkotas and Stergioulasl . Il999 



Lindblom. Mendell and Owen Ul999l ). The timescale for n = 1 polytrope from 



2 



gravity wave emission can be written as 
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where M, R, v represent the mass, radius and the rotation frequency of the 
star respectively. In general, the shear viscosity and surface rubbing dominate 
at low temperature, but in fact we can neglect the shear viscosity because 
the dissipation associated with a viscous boundary layer would greatly exceed 



2000; 
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(IBildsten and Ushomirsky 



The timescale due to surface rubbing is given 
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We now discuss the bulk viscosity coefficients of strange quark matter, hadron 
matter and their mixed phase matter to estimate the bulk viscous damping 
timescale in hybrid stars. The bulk viscosities are given approximately by 



(4) 



for strange quark matter (IMadsenl . Il992l ). and 



HP 



6 X 10" 



(5) 



for hadron matter (jlpser and Lindbloml . ll99ll ). where the constant k = 2/3, the 
angular velocity Q = 2ni> and p is the density, is strange quark mass. To 
construct hybrid stars, we consider the deconfiner nent transition occurring at 



high density ins ide the compact stars and use BP S (lBaym. Pethick and Sutherland 



197ll ) and GPS (lGhosh. Phatak and Sahul . ll995l ) equations of state for hadrons 
and MIT model for quarks. As in most cases, the variation of volume per unit 
mass and the chemical potential imbalance due to the oscillation in the star is 
so small that it is the liner part that contributes most to the viscosity, so we 



have used the relaxation time approxima t ion method(|Lindblom and Owen 



2OO2I : iDrago. Lavagno and Pagliara I . l2005l : iNayyar and Owen 1 . 12005| ) to com- 



pute the bulk viscosity of the mixed phase matter, 
^(7oo - 7o)^ 
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Fig. 1. The bulk viscosity of the mixed phase matter in the hybrid 



star for m. 



150Mev and 



moc 



del (|Ghosh. Phatak and SahJ . \l99^ ). 



165Mev in terms of GPS 



Here P is the pressure of the mixed phase, 700 and 70 are the "infinite" and 
"zero" frequency adiabatic index, u is the angular velocity of the perturbation 
and r is the relaxation time associated to the weak processes. As an exam- 
ple, we numerically calculate the bulk viscosity coefficient of the mixed phase 
matter which is showed in figure 1 for bag constant B^^'^ = 165 Mev. From 
numerical estimate, we easily find the fact that Cqp ~ Cmp ^ Chp- 

Now, we can use the bulk viscosities listed above to evaluate the bulk viscous 
damping times for the hybrid stars. For the given energy of modes E and its 
variation rate E = dE/dt, the timescale is defined as 
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The r-mo de energy in a rotating relativis t ic sta r with radius R has been ex- 
pressed asl lLindblom. Mendell and Owen I . Il999l ) 
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where a. is the normalization parameter, m is integer and p(r) is the density 
function ( Andersson and Kokkotas 1 . boOll ). The energy dissipation rate due to 
bulk viscosity can be estimated with the following formula 



bv 



C5a6a*d^x, 



(9) 



where 6a is expansion te nsor of Eulerian perturbation of velo city field, ex- 



pressed approximately as (lLindblom. Mendell and Owen I . Il999l ) 



S p 

6a = —iiui + mVl) — , 
P 



(10) 



here uj denotes angular frequency of the r-mode. In lowest-order approxima- 
tion, we easily find 5a ~ i?^(-^)'""'"^.In a stratiform star, we can decompose the 
energy E into four terms, i.e., the quark core, mixed phase, nuclei envelope 
and solid crust denoted respectively by QP, MP, HP and CR 

E = { j {Eqp + Emp + Ehp + Ecr) (11) 



with 



2 ^QP / \ 2m+2 

Eqp = ^{m + lf{2m + l)\Rt^pn' J pQp{r) [-^j dr, (12) 





o Rmp / \ 2m+2 

Emp = -^{m + l)'(2m + l)!Pjpfi' j pMp{r) ( j dr, (13) 

Rqp 

o Rhp / V 2m+2 

EHP = ^{m + lf{2m + l)\R%pn' j p^p(r) f -^j dr, (14) 

Rmp 



EcR='^{m + l)\2m+l)\R%pQ' j Pcnir) [^j dr. (15) 



Rhp 



A constant density in each phase, replaced by corresponding mean density, is 
a good approximation in terms of our calculations of the stellar structure(see 
figure 2). Considering Cqp,mp ^ Chp,cr, we have, according to (8) and (9), 

Eb. = {^Y'^-^Eqp + Emp), (16) 
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Fig. 2. The stratiform structure of the hybrid star with parameters M = 1.4M^, 
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9.87km, rus 



150Mev and B^/^ 



O' 

165Mev in terms of GPS 
model (|Ghosh. Phatak and Sahul . ll995 ). Here QP and MP represent the quark phase 
and mixed phase in the hybrid star respectively. 
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Thus, the timescale from Eq(4) is immediately expressed as 



(19) 



where A is a constant rela ted to the stellar structure that is showed in table 
1, T^^{R) can be given byf Andersson and Kokkotas 1 . 1200 ll ) 



(rik)^ (li^)"* (ix)"^ (is^)"* (T)^ 



Inserting (2), (3) and (18), we apply Eq(l) to evaluate the critical frequency of 
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Table 1 

The Values of A for different Structures of a Hybrid Star with M = IAMq 



51/4 [MeV] 


R[km] 


A 


165.0 


9.87 


0.6886 


175.0 


10.41 


0.6929 
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6.5488 
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Fig. 3. Comparison of r-i r iode instability windows of hybrid s tars w ith several other 
models (jAndersson et al.l . l200d : iMadsen I . l200d : Izheng et al.l . I2OO3I I. HS, SS and NS 
denote hybrid star, strange star and neutron star respectively, while the star symbol 
and the rectangle indicate the two observed fastest pulsars and the Low Mass X-rays 
Binaries. 

the rotating star as a function of temperature. Figure 3 shows the instability 
windows of the hybrid stars for M = I.4M0 and R ranging from 10~ 14km. 
As comparison, we also give the windows of strange star and normal star in 
figure 3. 

For hybrid stars, one notes that a pulsar will reach the frequency of about 
700 Hz( 1.43 ms) following a track coincident with the curve marking the 
instability region. So, the LMXB's are well within the region stable against 
r-mode instabilities, just spinning up the region unhindered by the instability 
due to accretion. For strange stars, there is the same mechanis m as hybr i d star 



model but a deviation of temperature from the pulsar data llMadsen I . 12000 
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Zheng et al.l . l2003l ). It is accepted that the neutron stars with crus t could 



also fit the pulsar data through r-mode runaway recycle mec hanism (ILevin 



1999 : Bildsten and Ushomirsky . 2000l : Andersson et al. . 2000 ). An important 
point is that, the fate of the crust of neutron stars is different from that of 
hybrid stars: the former is melted probably because of r-mode dissipation due 
to viscous boundary layer. We can imagine that the existence of the higher 
rotation pulsars is possible if there would exist larger quark matter core in 
the interior of a compact star. This may imply that submillisecond pulsars, if 
they exist, should be hybrid stars instead of neutron stars or strange stars. 

In conclusion, the bulk viscous damping timescale for hybrid stars is estimated 
for a range of parameters through which the r-mode instability window is 
obtained. The rapid rotation pulsar data is a significant indication of the 
existence of hybrid stars instead of neutron stars or strange stars. Naturally, 
this is a possible evidence of quark matter in the interior of the compact stars. 
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